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Appendix A 
Mark Ups for Claim Amendments 
1. (Amended) A method for predicting measurement data until a final 
time-point using given measurement data, En-w ^lcomprisin a th» * to pc 
of: 

a) "latching, using a processor , a stochastic process [is matched 
]to fthelsaid given measurement data; 

b) rxmjilncLsimulation runs of fthelsaid stochastic process^*© 
eamed^eut] from a giveri time-point until fth elsaid final time-point; 

c) re determining forecast measurement data[ is dotormincd] for 
each simulation run; and 

d) firedicting_measurement dataf-^red+eted] by stating a range 
of values, which is determined by [thejsajd forecast measurement data, and 
providing said predicted ma asurement Ha t a „ s useable n,,t p .,t 



m r 

■=«5- / 



2. /^mejidedLThe method as claimed in claim 1, [ifM^Jfurther 
comprising the ste ps re- 
determining a confidence range f> dotormincd ] for fth elsaid 
prediction of measurement dataf, whoro thol ; and 

e ' iminat ' n q a% lowest and b% highest forecast measurement dataf 
arc eliminated] . 

3. (Amended) The method as claimed in claim 2. fm lwherein 
[wWete]a% and b% are equal values . 

4. (Amended) The method as claimed in fono of thn preceding uUii„ b 
^Wefejftnejclaim 1, wherein said stochastic process is a non- 
homogeneous Poisson process. 

5. lAmendedjThe method as claimed in fono of thn preceding Ujj, IUj 

^efe]fthe]cl a im1, wherein said measurement data represents numbers 
of errors. 
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6. p^l f Amended) A method for predicting measurement data using 
given measurement data, [u^hichlcomprising the steps of: 

a) matching, using a processor, a stochastic process [ i o matched 
]to fthel said given measurement data;* 

b) [a range io asc e rtained, |^ orting [the-] probability values 
generated by [the]said stochastic process according to size, to provide a 
range around an expected value: and 

c) predicting measurement data[ io prodicted] within [the-]limits of 
fthelsaid ranr^ , and providing sai d predicted measurement data as 
useable output . 

! 7. (Amended) T he method as claimed in claim 6, [w-wh4eh] further 

comprising the steps of : 

[thel sorting said probability values generated by [the]said stochastic 
process[ aro ported] symmetrically by size around [the]sajd expected value. 

8. (Amended) A n arrangement for predicting measurement data until 
,.a,final time-point using given measurement data, [w^ereby-]comprisinq: 
' a processor ■■»» [ r rn wiHnri nn ^l. "^vina a CPU, bus, memory, and 
input/output controller, configured in such a way that: 

a) a stochastic process can be matched to [the]said given 

measurement data; 

b) simulation runs of the stochastic process can be carried out from 

a given time-point until the final time-point; 

C ) [the-]forecast measurement data can be determined for each 

simulation run: and 

d) measurement data is predicted by stating a range of values, 
which is determined by [the]said forecast measurement data,.?^ 
measurement data being output in a u seable form. 

; 9. (Amended) A n arrangement for predicting measurement data using 
given measurement data, [whnrobv Icomprising: 
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a processor unit [ i s providod and l . having a CPU, bus, memory, and 
input/output controller, configured in such a way that: 

a) a stochastic process can be matched to the given measurement 

data; 

b) a range can be ascertained by sorting probability values 
generated by [the]sajd stochastic process according to size around an 
expected value; and 

c) [thejsajd measurement data is predicted within the limits of the 
range, said measurement data being output in useable form 
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m ^ OV/ 868239 

W jWRec'dfgT/PTO 1 5jun 2001 

r D e scription l SPECIFICATION 

[ M e thod and arr a ng e m e nt for prod i ct i ng m e asuroment data using giv o n 

m e asurom e nt data ] 

TITLE 

METHOD AND ARRANGEMENT FOR PREDICTING MEASUREMENT DATA 
USING GIVEN MEASUREMENT DATA 

BACKGROUND OF THE INVENTION 

Field of the Invention 

FOOOn The invention relates to a method and arrangement for predicting 

measurement data using given measurement data. 

Description of the Related Art 

r00021 A technical system often requires facilities for forecasting based on 
known (measurement) data, particularly in the context of error susceptibility or cost 
estimates. 

r00031 Forecasts generated by experts are generally subject to errors. 

Experts cannot carry out exact analyses, at least of highly complex systems. 

r00041 A stochastic point process, in particular a Poisson process, is 
described in fWI Sidnev I. Resnick: "Adventures -ifv - in Stochastic Processes". 
Birkhauser Boston, 1992, ISBN 3-7643-3591-2, pp. 303-317 (Resnick) . 

Summary of the Invention 

rOOOSI The object of the invention is to allow the automatic prediction 
(forecast) of measurement data using given measurement data. 

r00061 This object is achieved in accordance with the [ features of tho 
ind e p e ndent patent claims. Deve l opm o nto l method and apparatus described 
below: developments of the invention are also described in the 
r d e p e ndont l followinq [etaimsl text . 

r00071 In order to achieve fthel this object, a method is f propos e d l provided 

for predicting measurement data using given measurement data, in which a 
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stochastic process is matched to the given measurement data. Simulation runs are 
carried out from a given time-point until a final time-point. The forecast 
measurement data is determined for each simulation run. Measurement data for the 
final time-point is predicted within a range of values, which is governed by the 
forecast measurement data. 

r0OO81 One development is to define a confidence range for the prediction of 

measurement data, where the a% lowest and b% highest forecast measurement 
data are eliminated. In particular, a% can equal b%. For example, a 95% confidence 
range can thus be defined by ignoring the 2.5% lowest and 2.5% highest forecast 
measurement data. 

[00091 One advantage is that the measurement data can be predicted 

(forecast) with an accuracy that is within a confidence range, from a given 
time-point. This makes it possible to identify- e^ . e.g.. the feasibility or 
impossibility of a task associated with the measurement data, at an early 
stage. Appropriate measures can therefore be initiated in order to counteract 
forecast impossibility. 

r0Q1Q1 [ On e advantage i s that tho measurement data can b o pr e dict e d 

(for e cast) with an accuracy that is within a confidonco range, from a g i ven t i me - 
point. Th i s makos it poocib l o to idont i fy ][-^gr][ tho feas i bil i ty or i mposs i bi li ty of a task 
associated w i th tho m e asur e ment data, at an early stage. Appropriat o measur e s can 
thoroforo bo i n i tiated i n ord e r to counteract forecast i mposs i b i lity. ]This is particularly 
important in the case of a complex system, e.g. A a software development process, 
where the extent to which a schedule can be followed before the software is 
completed can be shown in a subsequent test phase. Even more important in this 
context is the ability to adopt countermeasures at an early stage if a delay has been 
clearly identified, e.g^ in an integration test phase. This firstly affects the feasibility 
of the specified deadline (timescale) and secondly directly affects costs, since non- 
compliance with the agreed timescale often results in additional costs. 

room One refinement is for the stochastic process to be a non- 
homogeneous Poisson process. 

rooi21 In particular, the measurement data may in one refinement comprise 

numbers of errors. This applies to software development, for example, where the 
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level of maturity is documented in accordance with the errors measured in a test 
phase. Completion is directly dependent on this level of maturity. In other words, the 
software cannot be delivered to customers until most of the errors have been 
removed from the software. This is particularly important with regard to resources 
(required to test and correct errors) and costs (due to delayed delivery). 

£0013] In order to achieve the object of the invention , a method is also 
f propos e d l provided for predicting measurement data using given measurement 
data, in which a stochastic process is matched to the given measurement data. A 
range is ascertained, by sorting the probability values generated by the stochastic 
process according to size, around an expected value. Measurement data is 
predicted on the basis of this range, and in particular the probability values within the 
range. 

r0014l One development is for the probability values generated by the 
stochastic process to be sorted symmetrically by size around the expected value. In 
particular, this means that the highest probability value represents the middle of the 
range, i.e u the expected value, whereas the next highest probability value is 
arranged to the right or left of the expected value. The next highest probability value 
is then arranged symmetrically on the other side of the expected value, in turn. 

£00151 This analytical (design) procedure provides a range, where the breadth 
of the range in turn indicates which probability values are significant in the prediction 
of the measurement data. 

f00161 In one particular refinement, the breadth of the range is determined by 
ignoring the probability values that lie below a given threshold. 

r00171 This produces a range (confidence range), which has a specific 
breadth as a result of the threshold. This breadth corresponds to the certainty with 
which the measurement data is predicted. 

r00181 If one assumes that the stochastic process is a non-homogeneous 

Poisson process, then the non-homogeneous Poisson process defines a step size, 
particularly on a time axis t, which indicates when the next error will occur. One 
characteristic of the non-homogeneous Poisson process is that it has no memory, 
so that a "no-memory" search is carried out from each error that occurs at a specific 
time-point, for a time-point that indicates the next error. 
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[00191 In order to achieve the object of the invention , an arrangement is also 
r proposocfl provided for predicting measurement data using given measurement 
dataf , wh e r o bv l that has a processor unit [i s provid e d ]andjs configured in such a 
way that: 

a) a stochastic process can be matched to the given measurement data; 

b) simulation runs can be carried out from a given time-point until a final 
time-point; 

c) the forecast measurement data can be determined for each simulation 
run; and 

d) the prediction of measurement data for the final time-point can be 
predicted within a range of values, which is determined by the forecast 
measurement data. 

f00201 In order to achieve the object of the invention , an arrangement is 

further f propos e d l provided for predicting measurement data using given 
measurement data f, wh e r e bv l that has a processor unit ps prov i d e d land is 
configured in such a way that: 

a) a stochastic process can be matched to the given measurement data; 

b) a range can be ascertained by sorting probability values generated by 
the stochastic process according to size around an expected value; and 

c) the measurement data is predicted within the limits of the range. 

r0021l The arrangements are particularly suitable for carrying out the 
inventive method or the developments described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

f00221 Exemplary embodiments of the invention are shown and explained 
below with reference to the drawings, in which: 

Fig. 1 is a graph showing an accumulated number of errors over a test 

period; 

Fig. 2 is a graph showing the superimposed confidence ranges for different 

process models; 
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Fig. 3 is a flowchart showing the steps in a method for predicting 

measurement data using given measurement data; 

Fig. 4 is a further flowchart showing the steps in a method for predicting 

measurement data using given measurement data; and 

Fig. 5 fshewsl is a block diagram showing a processor unit. 

DETAILED DESCRIPTION OF THE INVENTION 

f00231 In order to be able to forecast a number of [e rrors to bo ]expected 

errors in a technical process, e.g. a in a software development process, non- 
homogeneous Poisson processes (NHPP) are calibrated[ T ] (i.e. A matched to 
measurement data, [er§r]such_as the occurrence of errors over time[ T ]) as follows: 

r00241 The following equation describes a counting process associated with 

the stochastic point process (non-homogeneous Poisson process): 

{N(t)} t€R + (1) 
r00251 n 

r00261 and a time-point to defines the end of a test period, i.e. x a time-point at 
which the given data ends. The stochastic processes 

{tf(t)} teR + and (2) 

{O(t)} R + (3) 

r00271 

f00281 are searched with 

p(u(t) < N(t) - N(t 0 ) * 0(t)| N(t 0 ) = n 0 ) £ a (4) , 

r00291 v v ' 

r00301 for all time-points where t>tO and given values a e. (0.1 ) (confidence 
level) and nO e. N. In particular, the following text examines the increases in the 
stochastic countings process in relation to the. time-point to. 

f003n In the present case, where equation (1) represents a non- 

homogeneous Poisson process, the following equation (cf. fWI Resnick ) 
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r00321 
r00331 applies for 



KKti) - »(t 0 ) = *) = «*(- Wi) - i(to)]) ■ fi(tl) " i(to)1 



0 £ to < ti < oo, £ e N 0 (6) 

T00341 

f00351 and an intensity (mean measure, mean value function) of 
r00361 

i: R + R + ,t f-» i(t>= EN(t) (7). 

r00371 Since the nature of the Poisson process dictates that the increases 
(error increases in this case) are independent of previous increases, equation (5) for 
the time-points t>to to define a (minimum) range 

[g U /g 0 ] s [guCtX 9o(t)] c n 0 <8) 

r00381 

1*00391 can be simplified to 

X p(N(t) - N(t 0 ) = /) > a (9) . 

roo4oi u 

f 00411 Due to the unimodal nature of the Poisson count density, a range [g u , 

g 0 ] can be determined as follows: 

f00421 Step 1 : Sort the elementary probabilities 



: = p(N(t) - N(t 0 ) = t\ £ e N 0 



r00431 

r00441 into descending order and label the values sorted thus using 



r00451 



P(0)' P(l)' • • • < i - e - {P0' Pl» • • •} = (p(0)' P(l)' ■ • •} 

p(o) * P(i) £•••>; 
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01 



10 



15 



20 



r00461 

roo47i 



Step 2: 
Step 3: 



^min : = min^ e N 0 

Determine 

Determine an index set 



Z P(i) * « 
i = 0 



f00481 

r00491 
rOQ501 

roosn 



where 



Pin' • • • ' Pi^ 



min 



} = |p(o) 



P (^in)J ; 



Step 4: Substitute g u := min{i} and h]g c := max{i} : 



The range from equation (8) is also referred to as the forecast range. 



Stochastic simulation (second approach) 

r00521 It is possible to determine the confidence range described using 

simulation, with the following steps: 

r00531 Step 1 : Start independent simulation runs based on the selected 
process model at time-point tO of the last error message m z N; 

r00541 Step 2: End a simulation run as soon as the required final time- 

point te is reached; 

r00551 Step 3: Repeat Step 2 until all simulation runs are finished; 

f00561 Step 4: Sort the numbers Ni(t e ) of the errors generated in the i- 

th simulation run in the time period (to, te), i=1, .... m, in descending order, and label 
the values sorted thus [pj^tj, . . . , N^fcJ : and 

r00571 Step 5: Substitute 



[00581 
[00591 



9u — fi (Lm-a/2j)( t «) and 
§ ° := S (rm-(l-a/2)l)( t e)' 



i.e. A eliminate the (100p]j_(1-a)/2)% lowest and highest values. 
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f00601 i ms produces the confidence range directly. 

r00611 Each individual simulation run is based on a simulation algorithm, 

which is known from (cf. ff^l Bratelv. et al.. 1987 ): 

10062] f Simu l atod l The simulated generation of intermediate arrival times for 
a non-homogeneous Poisson process is as follows : 

Stepl: Substituted : = sup {^(t)} , where: 

t>t- 



X(t) : = — 
dt 



(10) 



Step 2: Generate a (pseudo) random variable X that is exponentially 
distributed with the parameter X , Le. A x := - iog(u) / X , where U is equally 
distributed over (0, 1 )H; 

Step 3: Generate a random variable U that is equally distributed over (0,1 )H; 
and 

Step 4: If u < x(t s + x)/X, then substitute t* := t s + X; otherwise substitute t s 
:= t s + X and go to Step 1 . 

F00631 The example graph in Fig. 1 shows an accumulated number of errors 
during a given test period. From time-point to, it shows a prediction range Kl for all 
time-points to + x. 

r00641 The intensity i is normally derived from equation (10) for X . For 

example the result is as follows: 

r00651 a) ?i(t) = a . b • c . exp(- bt c ) . t 0 " 1 

r00661 (X(t) is strictly monotonously descending for [ , and un i modal for ] [ 
with a d e fin i tive max i mum at a po i nt ] c < i , and unimodal for c > i with a 
definitive maximum at a point 



r00681 b) Otherwise, X is derived in accordance with the above 
comments as follows: 
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f00691 

jj- = \*ks\ (c S 1) v (t s > t max ) 

[^(tmax) 

f00701 The graph in Fig. 2 shows the superimposed confidence ranges. In 
particular, this illustrates that possible forecasts become more scattered the further 
5 they extend into the future. In particular, confidence ranges calculated using 

different process models can be demonstrated in the same way as shown in Fig. 2. 

r00711 Fig. 3 shows a flowchart for the steps of a method for predicting 

^ measurement data using given measurement data. In Step 301 , a stochastic 

J| process, in particular a non-homogenous Poisson process (to represent a stochastic 

Co 10 count process), is matched to given measurement data. In Step 302, simulation runs 

are run from time-point to to a final time-point te that is to be forecast. In Step 303, 
=S for each simulation run, forecast measurement data is determined and a prediction 

O of measurement data is restricted to a range which is covered by the measurement 

g data determined by the simulation runs (see Step 304). In Step 305, a confidence 

15 range is determined in which a given proportion of the lowest and highest forecast 
ju measurement data is ignored in each case (this corresponds to the aforementioned 

range). The method terminates in Step 306. 

f00721 Fig. 4 shows a further flowchart for the steps of a method for 

predicting measurement data using given measurement data. In Step 401, a 

20 stochastic process, in particular a non-homogenous Poisson process, is matched to 
the given measurement data. Probability values are determined using the stochastic 
process, and these are sorted according to size around an expected value (see Step 
402). This sort operation results in the definition of a range, namely a confidence 
range in this case. The breadth of the confidence range is determined by comparing 

2 5 the accumulated probabilities with a given threshold. As described above, the 

confidence range gives a distribution or uncertainty, respectively, of a time-point to in 
the future, which allows the measurement data to be estimated in the future (see 
Step 403). The method terminates in Step 404. 
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100731 Fig. 5 shows a processor unit PRZE that mav be used to implement 
the inventive method. The processor unit PRZE comprises a processor CPU, a 
memory unit MEM, and an input/output interface IOS, which is used in different ways 
via an interface IFC: a graphics interface allows output to be viewed on a monitor 
MON and/or output to a printer PRT. Inputs are entered via a mouse MAS or a 
keyboard TAST. The processor unit PRZE also includes a data bus BUS, which 
provides the connection between a memory unit MEM, the processor CPU and the 
input/output interface IOS. It is also possible to connect additional components to 
the data bus BUS, e.g. additional memory, data storage (hard disk) or scanner. 

£00741 The C programming language is used in the following examples, which 
show an algorithm to define confidence ranges for forecasts and an algorithm for 
simulated definition of confidence ranges for forecasts. 

r00751 Program 1 : 

f00761 /* Definition of confidence ranges for forecasts */ 
I0077J I* based on the generalized Goel-Okomoto model 7 
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^include <stdiib.h> 
include <math.h> 
include <stdio.h> 



#define true 1 
#define false -1 



jj %3 
I = I 



double mv_genGO(doubl .double.double.double); double poisson(double,long); void ki nhpp(V 
int main(argc.argv) ~ v 

int argc; 
char *argvfl; 
{ 

double a,b,c.btst,kn; 
long low.upp.lauf; 

if (argc<7) { 

printf("\n\nZuwenig Argumente! \n\n M ); 

printffAufruf: %s <Par1> <Par2> <Par3> <Startzeit> <Endzeit>- 
"<KNiveau>\n\n\ argv[0]); return 1; 



O 

y2 a = atof(argv(1]); 

m b = atof(argv[2J); 

tfi c = atof(argv[3]); 

f» bt= atof(argv[4]); 

st= atof(argvI5]); 
kn= atof(argv[6]); 



for (lauf=1;lauf< ;lauf++) { 

kLnhpp(mv_genGO ( a,b t c.btbt+Iauf*(st-bt)/10. t kn f &low&upp)* 
pnntff 'Zeitpunkt: %8.2f Fehlerintervall: [%d,%d]\n- 
bt+lauf*(st-bt)/10., low, uppV 
) 

return 0; 

} 

double mv_genGO(x,a,b,c) 
double x,a,b,c; 

{ retum( a*(1.0-exp<-b*pow(x ( c))) ); } 

double poisson(lambda.wert) 
double lambda; 
long wert; 
{ 

long i; 

double itval.hv; 

if (lambda<600) { 
itval » exp(-lambda); 

for (i=wert;i>=i;H { itval •= lambda/(double)i; } 



r00781 
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r00791 



e!se{ 

hv= xp(-lambda/(double)weft); 
itval = 1.0; 

for (i=wert;i>=1;H { rtval •= lambda/(doub! )i*hv; } 

return ( rtval ); 
} 

void kLnhpp(mv_nhpp, parl_nhpp, par2_nhpp, par3_nhpp, 

startjime, stop_time, k_niveau t lower, upper) double mv_nhpp(doubIe.double.double.double); 
double par1_nhpp, par2_nhpp, par3_nhpp, start.time, stop.time, k_niveau: long lower, 'upper 

long lauf; 

int Iborder,modJow,mod_upp; 
double sum,tmp_mv, val_l ( val_u; 

tmp_mv = mv_nhpp(stop_time t par1_nhpp t par2_nhpp,par3_nhpp) - 
mv^nhpptstart^me^arl^nhpp.pa^nhpp.parS^hpp); lauf = (long)tmp_mv; 
*lower ~ lauf-1 ; 

•upper = lauf+1; modJow= false; mod_upp= false; sum = poisson(tmp_mv,lauf); valj = 

poisson(tmp_mv t *lower); val_u = poisson(tmp_mv,*upper); 
while (sum<k_niveau) { 
if (valj<val_u) { 
sum += val_u; 
(*upper)++; 
(border = false; 
mod_upp = true; 

val_u = poisson(tmp_mv,*upper); 
} 

else { 

sum += vaM; 
Power)-; 
Iborder = true; 
mod_low = true; 

valj = poisson(tmp_mv,*!ower); 
} 

} 

if (Iborder == true) { (*lower)++; } 
else { Cupper)-; } 

if (modjow = false) { (*lower)++; } 
rf (mod_upp == false) { (*upper>-; } 

return; 
} 

Program 2: 



r00801 



/* Simulated definition of confidence ranges for forecasts 7 



roo8ii 



/* based on the generalized Goel-Okomoto model */ 
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#»ndude <stdlib.h> 
include <math.h> 
#include <time.h> 
#indude <stdio.h> 
#include <values.h> 

#define true 1 
#define false -1 

double drand48(void); 
void srand48(long); 

double sim_exp(double); double lambda_genGO(double,double,double,double); void sim nhpp(V 

int main(argc.argv) 

int argc; 

char *argvQ; 

{ 

time_tt; double a,b,c,btst ( pnt[1000000] t check_time{12J; long lauf.no ontseed run; intclaut 
FILE *datei; 
if (argc<6) { 

printf("\n\nZuwenig Argumente! \n\n'); 

printfCAufruf %s <Par1> <Par2> <Par3> <Startzeit> <Endzeit>\n\n - , 
argv[0]); return 1; 
} 

datei = fopenrsim.seedYr"); 
if (datei==NULL) { 
seed_njn = 1 ; 
} 

else{ 

fscanf(datei t "%6d" ( &seed_run); fclose(datei); seed_run++; 



datei = fopen( ,, sim.seed^"w+ ,, ); 
fprintf(datei. "%6d\n\ seed run ); 
fclose(datei); 

time (&t) ; r Initialisierung des */ 

t += seed_run*100 ; /* Zufallszahlengenerators 7 

srand48 ((unsigned long) t) ; r mit Hilfe der Systemzeit 7 

a = atof(argv[1]); 

b = atof(argv[2]); 

c = atof(argv[3]); 

bt= atof(argv[4]); 

st= atof(argv[5]); 

sim_nhpp(lambda_genGO,a,b,c,btst&pnt ( &no_pnt); 
for (lauf=l;lauf«no _pnt;iaufw-) { 
printff%l5.7f %lOd \n\ pntjlaufl, lauf); 
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f00831 



datei = fbpenfki.tmpYa B ); 
for (Iaufs1;lauf< ;taufM-) { 
check^timeflaLrf] = bt+lauT(st-bt)/10.; 

check_time[1 1 J = pnt{no_j>nthM; r groOer als die grOBte 
simuiierte Zeit •/ 
ciauf = 1; 

for (lauf=1;lauf<=no_pnt;lauf++) { 

while (pntIlaufJ>=check_tiniefciaufD { fprintQdatet. "%8.2f %6d check_time[clauf], lauM); clauf~; 
} 

if (pntfno_pntJ <check_time[10]) { 

for {lairf=dainiauf< : ia~uf*+) { fprintf(datei, "%8.2f %6d check_tirne[lauf]. nojnt); 
} 

fprintf(datei, "\n'); 
fciose(datei); 

return 0; 
} 

double sim_exp(lambda) 
double lambda; 

{ retum( -log(drand48<))/lambda ); } 

double lambda_genGO(x,a 1 b,c) 
double x,a,b,c; 

{ return( a*b*c*pow(x.c-1 )*exp(-b*pow(x ( c)) ); } 

void sjm_nhpp(lambda_nhpp, part_nhpp, par2nhpp, par3 nhpp 

^ir^nhn^n^ 6, no ^°^ doub,e ^mbda.nhp^double.double.double.double); double 
^ pan_nhpp, par2_nhpp, par3_nhpp, start.time, stop.time; double pathQ; long -no joints; 

double simjime,x,u,x_bar,lambda bar 

*no_points=0; 

srm_time = startjime; 

do{ 

if <par3_nhpp<=i) { ,ambda.bar * lambda .nhpp(sim_time,par1_nhpp,par2,nhpp,par3_nhpp); 
else { 

x bar = Pow((par3.nhpp-1.0)/par2.nhpp/par3_nhpp t l.0/par3 nhppV 
» (sim_time>=x_bar) { ~ 

lambda.bar = Iambda_nhpp(sim_time.parl.nhpp.par2_nhpp,par3.nhpp); 

else { 

lambda.bar = lambda.nhppCx^bar.parl.nhpp.pa^nhpp.parS.nhpp); 
} 

x = sim_exp(lambda bar); 
u = drand48(); 



if (" <= lambda.nhpp(sim_time+x f par1.nhpp ( par2_nhpp f par3jhpp^ambda bar) { fno _pointsH 
path['no_points]=sim time+x; ~ 

} 

srm_time+=x; 
} 

while (sim_time<=stop_time); 

return; ~" 

) 
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Program 3: 

/* Definition of confidence ranges from the simulation data */ 
/* (the simulation data is sorted into ascending order) */ 
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include <stdlib.h> 
#incJud <math.h> 
#indud <stdio.h> 



int qsortJcmp(int*,int*); 
int qsortjcmp(x.y) 
int *x, *y; 
{ 

if (*x<*y) { return ( -1 ); } 
else if (*x—*y) { return ( 0);} 
else { return ( 1 }; } 

int main(argc,argv) 
int argc; 
char *argvQ; 
{ 

int pnt(11][1 00000]; 

int qs[100000J; 
— char *dname; 

"~ int fracj; 

^ long lauf ( lower_bound,upper_bound; 

EQ long l,no_pnt,seed_run; 

ffi double ctime[1 1 J.x" 

Ol FILE *datei; 



if (argc<3) { 

prin *i^ W): PrintfTAuftuf. %s <Dateiname> <Konfidenzniveau (in 

^ /o /o)>\n\n p argvfOJ); return 1 ; 



dname = argv[1J; 

frac » 100-atoi(argv[2]); 

lauf=0; 

datei = fopen(dname, V); 
if (datei==NULL) { return 1; } 
else{ 

while (!feof(datei)) { 
lauf++; 

for (i=1;i<=9;i++){ 

fscanf(datei,-%8lf%6d". &ctimep], &pnt[ij[lauf]); 
fscanf(dater/%8lf %6d \n", &ctime[10], &pnt(10][lauf]); 

fclose(datei); 
} 



lf°)K 0 K nd = (^ n ?)floor{laurfrac^200.); upper.bound = (long)ceil(lauf(200..frac)/200 )• 
rf (lower_bound<1){lower_bound= 1;} ' W,A 

f00871 
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print^\jio/ o2d o /o o /o . Sichertleitsbereich ^ %d S imulationsiaufen\n\n- 
1 OO-frac.lauf); 
for (i=.l;i< 

for (»=1;l<=lauf;f++){ 
qs[Q = pnt!iJ[IJ; 
> 

qsort(&qs[1], lauf, sizeof(int), &qsort_icmp); 
pnntffZeitpunkt: %8.2f Fehlerintervall: t%d,%d]\n" 
ctimefi], qs[lower_boundJ, qs[upper_boundJ); 



return 0 ; 
> 



mm The above-described metho d and apparatus are illustrative oftho 

principles of the prese nt invention. [Rofcronoos:] 

[f4f 



-][ Sidnoy I. Rosnick: "Adv o nturoG ][4ft-][S t o chastio PrococooG", 
Birkhaus o r Bocton, 1002, ISBN 3 7613 35 01 2, pp. 0 0 3 3 17 ] 
[{29 Rrato l y ot a l ., 10 0 7] 

j!!".!!!! r0 "t m0 ^ fiCat J 0nS 3nd adaptations wi " be readily ap pa rent to the,** 
skilled in this art withoutdepartina from t h e spirit and . ^ppe of the p r^nt 

invention 
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[Abstfaet] 

[M ot h n d and arr a iNj , - mont f o r p rou.bU. hj i t . . wu , 0 ,„ u 
m on nnromont UdLi l ABSTRACT 

I2Q891 A method and an arrangement are [pfepesed! provided for predicting 
measurement data using given measurement data, in which a stochastic process is 
matched to the given measurement data. Simulation runs are carried out from a 
g.ven tone-point until a final time-point. The forecast measurement data is 
determined for each simulation run. Measurement data for the final time-point is 
predated within a range of values, which is determined by the forecast 
measurement data. 
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